Diffuse emission in the presence of an inhomogeneous spin-orbit interaction for the 

purpose of spin filtration 
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A lateral interface connecting two regions with different strengths of the Bychkov-Rashba spin- 
orbit interaction can be used as a spin polarizer of electrons in two dimensional semiconductor 
heterostructures [Khodas et ai, Phys. Rev. Lett. 92, 086602 (2004)]. In this paper we consider 
the case when one of the two regions is ballistic, while the other one is diffusive. We generalize the 
technique developed for the solution of the problem of the diffuse emission to the case of the spin- 
dependent scattering at the interface and determine the distribution of electrons emitted from the 
diffusive region. It is shown that the diffuse emission is an effective way to get electrons propagating 
at small angles to the interface that are most appropriate for the spin filtration and a subsequent 
spin manipulation. Finally, a scheme is proposed of a spin filter device, see Fig. [51 that creates two 
almost fully spin-polarized beams of electrons. 
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PACS numbers: 72.25.Dc, 72.25.-b, 72.10.Bg, 73.63.-b 



I. INTRODUCTION 

Recently we proposed^ to use a lateral interface be- 
tween two regions of the two-dimensional (2D) electron 
gas with different strengths of the Bychkov-Rashba 2 spin- 
orbit (SO) interaction as a spin-polarizing element for 
the purposes of spintronics^*^ The lateral interface in- 
troduces the space-varying SO interaction which leads 
to spin-dependent refraction of spin carriers passing 
through the interface. 5 Consequently, an electron beam 
with a nonzero angle of incidence after passing through 
the interface splits into beams with different spin polar- 
izations propagating in different directions (see Fig. ^ for 
the scattering at the lateral interface between the regions 
with and without SO interaction). Using an interface 
with an inhomogeneous SO interaction as a principal el- 
ement of spin-based devices we outlined the schemes for 
a spin filter, spin guide, and a spin current switch (spin 
transistor). This program promises to build spintronics 
devices avoiding magnetic materials that are not conven- 
tional for the semiconductor industry. It can be realized 
in the gated heterostructures with a sufficiently strong 
Bychkov-Rashba SO interaction^**^ by manipulating the 
gates i 

The effect of the separation of the trajectories of elec- 
trons with different spin polarizations (chiralities) influ- 
enced by the space-varying SO interaction has the same 
grounds as the double refraction (birefringence) of light 
in uniaxial crystals exploited in optical devices for the 
polarization of light. The separation of the trajectories 
of carriers of different spin has been observed recently*^ in 
the case of a homogeneous SO interaction as a result of 
action of a perpendicular magnetic field. The separation 
of the trajectories after reflection at a lateral potential 
barrier in the presence of the SO interaction has been 
observed in Ref. 

Two facts that can be useful for the purposes of the 
spintronics have been found in Ref. U in the analysis of 
the spin-dependent scattering of electrons incident on a 



lateral interface with a SO interaction varying in the di- 
rection normal to the interface. First, there exists an 
interval of outgoing angles 8 C < 6 < it/2 where only 
electrons with a definite spin chirality can penetrate (see 
Fig. . If it is possible to collect electrons from this in- 
terval, one will have an ideal spin filter. Second, electrons 
of this chirality exhibit a total internal reflection for an 
angle of incidence in an interval tp c < tp < tt/2, where ip c 
is a critical angle of the total internal reflection. It is clear 
from these observations that the electrons propagating at 
small angles to the lateral interface are most sensitive to 
the variation of the magnitude of the spin-orbit interac- 
tion, and therefore such electrons are most appropriate 
for spin control and manipulation. Hence one has to find 
a way to create (and collect) flows of electrons of high 
intensity that are almost tangential to the interface. In 
this paper we suggest using the diffuse emissionii*^ as a 
possible solution of this task. 

We show that making one of the two regions connected 
by the lateral interface to be diffusive is an effective way 
to achieve a flat angular distribution of particles emitted 
from the diffusive region into the clean one. The effect of 
flattening of the angular distribution of the emitted elec- 
trons is the robust property of the diffuse emission which 
holds despite the presence of a spin-dependent reflection 
at the interface. Due to the flatness of the angular dis- 
tribution a substantial portion of electrons propagates 
at small angles to the interface and these electrons are 
suitable for the spin filtration and subsequent spin ma- 
nipulation. In Fig. a scheme of a device is presented 
which operates as a spin filter with a high level of spin- 
polarization of a filtered current. 

The paper is organized as follows. In Section [H] we 
present the resultsi of the analysis of scattering of elec- 
trons at the lateral interface between the two regions with 
different magnitudes of the Bychkov-Rashba term. In 
Section ITTT1 we analyze the transport of the 2D electrons 
in the presence of the Coulomb interaction near the in- 
terface between the diffusive and ballistic regions. In 
Section lTvl we reconsider the problem of the diffuse emis- 



2 



sion (Milne problemAiii£) for an arbitrary dimension. In 
Section [V] we generalize the technique developed for the 
problem of the diffuse emission to the case of the spin- 
dependent scattering at the lateral interface. We use a 
semiclassical approach with electrons moving along the 
classically allowed trajectories when the spin of electrons 
is the only element treated quantum-mechanically. For 
this purpose, we analyze the transport near the interface 
in terms of the spin density matrix. In the Summary 
a scheme of a spin filter device that creates two beams 
of electrons of a very high level of spin polarization is 
presented. 



II. SPIN-DEPENDENT SCATTERING AT THE 
LATERAL INTERFACE 

Consider a two-dimensional electron gas confined in 
the xz plane by a potential well in the semiconductor het- 
erostructure. Generally, the potential well has the shape 
of an asymmetric triangle, and, consequently, there is 
a direction of asymmetry, 1 , perpendicular to the elec- 
tron gas plane. This leads to the appearance of the 
specific spin-orbit interaction terrn^ in the Hamiltonian, 
a(p x l)er. We will consider the case when the parame- 
ter a varies along the x direction. The direction of 1 is 
chosen as 1 = — y. Generally, the Hamiltonian has the 
form: 



Hr 



-(1 x cr)[a(x)p + pa(x)]. 



(1) 



Here B(x) describes the varying bottom of the conduc- 
tion band which may be controlled by gates. The current 
operator corresponding to this Hamiltonian contains a 
spin term: J = ^ + a(x)(l x <x). The presence of spin in 
the current operator implies that in the process of scat- 
tering at the lateral interface with varying a the conti- 
nuity conditions for the wave function involves the spin 
degrees of freedom of the electrons. This makes the elec- 
tron scattering at the interface to be spin-dependent. 

To diagonalize the Hamiltonian with the Bychkov- 
Rashba term in the regions of constant a one has to 
choose the axis of the spin quantization along the di- 
rection (1 x p). Then for the two chiralities (referred to 
below as " + " and " — " modes) the dispersion relations 
are given by 



E ± = — ± ap 
2m 



B. 



(2) 



Correspondingly, the momenta of the waves of a given 
energy E involved in the scattering at the interface be- 
tween the two regions with different a are determined as 
follows: 



Pso = m(^2(E - B)/m + a 2 Ta) 
= mvp(V 1 + a 2 T 5). 



Here vf = \/2(E — B)/m, and we introduce a small 
dimensionless parameter a = a/vF which we will use 
throughout the paper. Notice that at a given energy 
the velocity of electrons of both chiralities is the same: 

v — dE^ I dp = vfv1 + a 2 . 

Let us discuss the kinematical aspects of the electron 
scattering at the lateral interface. An incident (nonpo- 
larized) beam comes at angle ip from the region without 
a spin-orbit term (or with a suppressed spin-orbit term) , 
denoted as the N region, and when transmitted into the 
SO region splits into two beams of different chirality that 
propagate at different angles 6' ± . Figure illustrates the 
scattering for the simplest case when a(x < 0) = 0. The 
conservation of the projection of the momentum on the 
interface together with Eq. J3J) determine the angles of 
the transmitted and reflected beams (SnelPs law): 



Pn sin ip 



p^O sin t 



(4) 



where pn is a momentum of an electron in the N region 
and pg are the momenta in the SO region after passing 
through the interface. From Eqs. © and (0J it follows 
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FIG. 1: (Color online) The refraction of electrons at the inter- 
face between the regions with (SO) and without (N) spin-orbit 
interaction. A beam incident at angle ip splits after the re- 
fraction into two beams with " + " and "— " chiralities that 
propagate at angles 8 ± (denoted by red and blue colors, re- 
spectively), (a) 9 C determines the limited aperture for " — " 
chirality; in the angle interval 8 C < 8 < tv/2 only electrons of 
+ chirality can penetrate, (b) <p c is an angle for total internal 
reflection for electrons of the + chirality. 

that the SO region is more "optically" dense for the + 
mode (i.e., it has a smaller wave vector) and less dense for 
the — mode. Correspondingly, the + mode is refracted to 
larger angles than the — one, and therefore the + mode 
exhibits a total internal reflection with a critical angle 
tp c . As to the — mode it has a limited aperture in the 
SO region for outgoing angles: 9 < C < n/2. 

Remarkably, the angle interval where only the + mode 
can penetrate is not narrow as its width has a square root 
dependence on 5. It follows from SnelPs law that (tt/2 — 
8 C ) w (■k/2~<p> c ) V2a. Actually, one can reduce 6 C even 
further. With the gates acting selectively on the different 
regions of the electron gas, SB = B(— oo) — B(+oo) ^ 0, 
one can alter the position of the bands relative to the 



3 



Fermi level in the N and SO regions. A simple analysis 
shows that with an increase of 5B (i.e., lowering pp in 
the normal region) the angle interval (tt/2 — 9 C ) grows 
and reaches 1\fa. Starting from this moment the angle 
interval suitable for spin filtration narrows and eventually 
becomes ~ 5, instead of ~ \fot. 

The problem of scattering of electrons at a lateral in- 
terface between the two regions with different magni- 
tudes of the Bychkov-Rashba term has been considered 



in Ref. Q for the two limiting cases of a sharp, X/d > 1, 
and smooth, X/d <C 1, interface, where A is an electron 
wavelength and d is an effective width of the interface. 
Here a qualitative description of the analysis of the spin- 
dependent scattering at the interface will be presented 
only. 

A scattering state of an electron coming from the N 
region in the state e % u>* x +P* z )x+ ' IS given by 



where X^/so are s P mors corresponding to the ± chiral 
modes in the N and SO regions, and r and t are the 
amplitudes of the reflected and the transmitted waves. 
In Eq. J2J the interface is at x — and for simplicity we 
limit ourselves to the case of the interface between the N 
and SO regions. A similar expression holds also for 
which evolves from the incident state Xn • 

The flux of particles impinging on the interface at a 
given angle (p can be defined as I^dedtpds, where ds is 
the cross-sectional area (width) of a beam. The inten- 
sity lf°of the transmitted flux of a given chirality can 
be found from the relation \Tso<— n(<P)\ I ^ cos ipdip — 

jSO 

cos 6 d9. Here the cosine factors take into consid- 
eration the change of the cross-sectional width of the 
beam as a result of the scattering, while \T'so->—n( 1 p)\ = 
\t\ 2 (jij /i)f) (we use the fact that for the two chiralities 
the velocities are the same). Finally, the intensity of the 
refracted flux relative to the incident one is equal to 

i!°/if = {de/d v )-\v so /v N )\t\\ (6) 

For a sharp interface the amplitudes r, t can be found 
from the continuity conditions, while for a smooth inter- 
face one can conduct the analysis of the refraction using 
a small parameter r] = (da/dx)/apF ~ X/d <C 1. In the 
latter case the electron spin adjusts itself adiabatically to 
the momentum keeping its polarization in the direction 
perpendicular to the momentum. 

It has been shown in Ref. Q that in the course of re- 
fraction at the interface with a <C 1 transitions between 
states of different chirality are strongly suppressed. For 
that reason, the drop of the intensities of the transmitted 

electrons without change of their chirality t ++ and t , 

see Fig. |21 , occurs practically only due to the reflection 
which becomes decisive only for (p > ip c . In particular, for 
a smooth interface the probability of the reflection out- 
side the region of the total internal reflection is almost en- 
tirely suppressed because the matrix elements describing 
reflection are integrals of the rapidly oscillating functions. 
Consequently, the transmission amplitudes presented in 
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7T/4 <p c 7T/2 



FIG. 2: (Color online) The intensities per unit outgoing angle 
of the electrons transmitted without change of their chirality 

~ (dO + /dtp)~ 1 \t++ | 2 and~ {dd~ /dip)' 1 ^ | 2 per unit angle 

as a function of an angle of incidence for sharp (solid line) and 
smooth (dashed line) interfaces (Ref. 

Fig. have almost rectangular shape for a smooth inter- 
face. 

Summarizing the results obtained in Ref. |jj for the 
cases of the sharp and smooth interface one can state 
that for both the discussed cases an electron in a state 
of a definite chirality propagates along the classically al- 
lowed trajectory for this chirality, while a change of the 
chirality is very ineffective. 

The case of the N-SO interface described so far was 
taken mostly for illustration. Actually, any interface with 
a change of a results in splitting of the trajectories that 
can be used for the purposes of spin polarization and 
filtration. In the analysis above, one should replace a by 
the difference of the strength of the spin-orbit interaction 
Saso across the interface. 

III. KINETIC EQUATION, 
ELECTRONEUTRALITY, AND BOUNDARY 
CONDITIONS 

In this section we discuss the kinetics of the two- 
dimensional (2D) electron gas near the boundary between 
the diffusive and ballistic regions. We show that within 
the linear response approximation the equation determin- 
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ing the current distribution function and the equation 
determining the density and potential profiles are decou- 
pled, except for boundary conditions. This analysis justi- 
fies the concept introduced by Landauer*i4 to determine 
the transport properties it is enough to consider noninter- 
acting quasiparticles that propagate without interaction, 
i.e., ignoring any effects related to the redistribution of 
the potential due to the Coulomb interaction of electrons. 
The problem requires, however, fixing up the boundary 
conditions, which we perform below for a particular ge- 
ometry. 

The kinetic equation describing a stationary flow of 
electrons by a distribution function n p (x) is 



de p (r) dn p (r) _ de p (r) dn p (r 
dp dr dr dp 



St{n p (r)}. (7) 



The electron flow is forced by the electric field, 
de p (r)/dr = eV$(r). Correspondingly, Eq. J7J should 
be supplied by the Poisson equation which when limited 
to the 2D plane is 



A$ = -4Trep 2D S{y). 



(8) 



Here the Laplacian acts in 3D space, while p 2D is the 
deviation of density of the 2D electron gas from the equi- 
librium value. Therefore Eq. (jSJl should be supplemented 
with the 3D conditions on <f> that will take care of the 
3D environment of the 2D electron gas. This introduces 
an element of nonuniversality into this problem. For- 
tunately, in the linear response approximation the sit- 
uation is much more tractable and the problem of the 
current flow of the 2D gas becomes self-contained (see, 
e.g., Ref.[H). 

In the linear response approximation one keeps terms 
linear in V$(r) and 8n p (r) only: 



d5n p (r) 
dr 



eV$ r v F — £ 

de 



St{Sn}, (9) 



where 8n p (r) = n p (r) — wP F and n° F is the Fermi-Dirac 
equilibrium distribution. (We assume throughout the pa- 
per that the spatial variation of the Fermi-energy level as 
well as other parameters of the electron liquid such as the 
density of states and the screening length occurs on dis- 
tances greatly exceeding the wavelength and the mean 
free path. For that reason, the spatial variation of the 
parameters characterizing the electron liquid is ignored 
below.) In the presence of the electron-electron interac- 
tion there is an additional force that originates from the 
interaction of the quasiparticles. This effect has been ac- 
counted for by a substitution of the local equilibrium dis- 
tribution Sn p (r) in place of Sn p (r) in Eq. ©, see Ref.Url 
Following Refs. D]J and E one can shift the distri- 
bution function by a local value of the potential $(r) 
introducing a displacement function f(r, ip>) 



Sn p (r) 



/(r,¥>)], 



(10) 



where the direction of the momentum p is given by the 
angle (p. In terms of the displacement function /(r, <p) 
the system of Eqs. Q and acquires the form 



dr 



= St{/}, 



(11) 



A$ - 2n 2D ^(r)8(y) = -2k 2D </(r, <p)) v 8(y), (12) 

where k 2 d — 2ire 2 dn/dp is the inverse screening 
length of the 2D electron gas, which in the Thomas- 
Fermi approximation is equal to 2e 2 m; (f(r,ip))^ = 
f (dip / 2ir) f (r , ip) is f(r,ip) averaged over directions of 
momentum. Thus, although Eq. together with the 
Poisson equation JSJ constitutes a system of two coupled 
equations, the potential <£>(r) drops out from Eq. i|ll|) 
governing the function f(r,ip). Given boundary condi- 
tions, the displacement function /(r, ip) determines the 
current distribution without any feedback from Eq. 112|) : 



j(r) = 2e 



d 2 p 



v F 5n p (r) = 2e 2 v 2D (v F (ip)f(r,tp)) 



v> ' 

where v 2F) is the density of states per one spin species, 
V2D = m* /2tt. 

When one is interested in the connection of the current 
density with the distribution of the potential and density 
in the 2D electron gas, Eq. Ijl2|l should to be involved 
along with the relation connecting p 2D with $(r) and 



d 2 p 



Sn p (r) 



(2n) 2 

= (KW27re)[</(r,p)>. 



$(r)]. (14) 



For good enough conductors the Poisson equation re- 
duces to the condition of the electroneutrality 18 (which 
is valid in any dimension). Under these circumstances 
typical distances on which the potential $ in the left- 
hand side of Eq. I|12l) changes are much longer than 
the screening length k^ d . Correspondingly, |V$| <C 
| K 2L>(/(r, ^)) |) | K 2D ( I > (r)|, and the Poisson equation 
reduces to the condition of the electroneutrality: 



<Wr)=0; (f(r,<p)) = #(r). 



(15) 



Now we turn to the discussion of the boundary con- 
ditions. Let us assume that the diffusive region has a 
stripe geometry with the x axis directed along the stripe 
(see Fig. |3J). In the diffusive region the collision term 
is controlled by the elastic relaxation time r e i, and the 
kinetic equation acquires the form 



dfjr, p) 

VF ^ = &t{ J j elastic 



dr 



/(r,y)-(/(r,y))„ 

Tel 



(16) 



Here we assume that the impurity scattering is of short 
range nature. (The angular dependence of scattering of 
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electrons by the impurities does not influence the distri- 
bution function of electrons deep inside the diffusive re- 
gion, but is important for the angular profile of the diffuse 
emission.) With the elastic mean free path I = vpT e i used 
as a unit of length, the kinetic equation can be rewritten 
in terms of the dimcnsionless variables £ = x/l: 



- C0Sl p—f + = (f(r,<p)) v . (17) 

Here the minus sign in the first term containing d/dC, 
appears because ip is chosen as an angle formed by a 
momentum with the direction —x (note that —x directs 
inwards the diffusive region). 

In the current carrying state electrons deep inside the 
diffusive region are distributed according to the Drudc 
form. Correspondingly, at a distance (counted from the 
interface) exceeding few I Eq. (|17|1 has a solution: 

/Dr(C^) = -^ 1 (cOS^ + C). (18) 

Here the factors are fixed in such a way that the current 
J(C) defined as 

J(0=-2tt (cos v>/(C, (19) 

is equal to J for the distribution function /Dr(C)V)- 
With the use of Eq. I|13|) the physical current j = j x = 
—2e 2 V2DVp (cos</?/(C, </?)) can be related to the current 
J defined in Eq. (|T§)| as follows: 

j = PF (e 2 /2rr 2 )J. (20) 

To obtain the relation between the current j and the 
electric field let us apply a gradient to the both sides 
of Eq. (|14fl with the distribution function /b r used for 
f(r,p). Then, together with the fact that the electric 
field E = — V<&, one gets the relation: 

j = oE- D ch V(eSp 2D ), (21) 

where the conductivity a — 2e 2 ^2D(^ 2 /2r) while the 
diffusion coefficient of charge density, D c h, corresponds 
to the Einstein relation, D cn = (dp/ dn)a / e 2 . Under 
the condition of the electroneutrality Sp 2 r>{r) = an d 
j = aE. [Actually Eq. H21|l does not require the Drudc 
form, /□,., of the distribution. It holds in the diffusion 
approximation when the distribution function has only 
two first harmonics f(x, p) — p Q >(x) — f^(x) cosp and 
the functions f^°'(x) and f^(x) vary on scales much 
exceeding the mean free path I.] 

To proceed further with Eq. (|llf> . one has to specify 
the distribution of particles incident from the terminals 
located on the ballistic side of the device under discus- 
sion. Generally, this distribution is not universal as it 
depends on a particular geometry. For definiteness, we 
consider the case when the diffusive region that has stripe 
geometry runs into a ballistic basin, see Fig. [3| (a stripe 




FIG. 3: The physical setup which leads to Eq. (12 H as a 
boundary condition for the kinetic equation. The diffusive 
region has a geometry of a stripe that runs into a ballistic 
basin. The current density vanishes inside the ballistic region 
at a typical distance L <C U n . Trajectories that start from 
a terminal on the ballistic side of the device and reach the 
diffusive stripe are shown by dashed lines. 

is wide enough allowing the transverse quantization of 
electrons to be ignored) . Following Ref . [la, we use the 
method of characteristics to determine a distribution of 
particles that impinge onto the diffusive stripe from the 
ballistic region (corresponding family of angles ip will be 
denoted as p b ^ d ). 

In the ballistic region the collision term is controlled by 
the inelastic relaxation time Ti n , and the kinetic equation 
acquires the form 

df(r,p) _ 

V F — f inelastic 

= _/('.*>)-*('), (22 ) 

Tin 

The solution of Eq. (|22p for the directions (f b ^ d is given 
by 

f(T, V )= / dt *F(i)r, V ], (23) 

J~ oo T in 

where r(t) TtV is the trajectory (the characteristics) of an 
electron that starts at the remote past from a terminal 
on the ballistic side of the device and reaches the point 
r with the momentum directed along ip at the moment 
t = 0; see dashed lines in Fig. We assume that the cur- 
rent carrying area widens sharply on the ballistic side of 
the setup. Correspondingly, the current density vanishes 
inside the ballistic region at a typical distance L. Under 
the condition = VFi~i n , the integral in Eq. H23JI 

is accumulated at distances where electrons are at the 
equilibrium and the potential <&(r) is equal to its equilib- 
rium value at a terminal deep inside the ballistic region, 
$(+00). Then, it follows from Eq. (|23|l that for incoming 
directions f(r,p b ^ d ) = <f>(+oo). This provides us with 
the boundary condition to be imposed at the interface, 
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0, on the incoming part of the function: 

/(X = 0,^^) = *(+«>)• 



(24) 



The obtained boundary condition is isotropic. This re- 
markable feature is a consequence of the choice of the 
proper geometry. 

Together with the current distribution deep in the dif- 
fusive region given by Eq. (JTSJ, the relation l|24|) consti- 
tutes the full set of the boundary conditions needed for 
the solution of Eq. I|16(l . Two remarks are now in order 
to complete the discussion of the boundary conditions. 

(i) Notice that $(x = 0) ^ $(+oo). The point is 
that due to the abrupt change in scattering the solu- 
tion /(r, (f) of Eq. H16(l has a singular derivative near the 
interface. Under these circumstances, one cannot neglect 
the term A$ in Eq. 1|12|) . and therefore the condition of 
the electroneutrality is violated in the vicinity of the in- 
terface in a strip of a width oc /t2£> • The deviation of 
the density distribution Sp 2 o ( r ) from the equilibrium at 
the interface leads to the variation of the potential, and 
hence &(x = 0) differs from <&(+oo). 

(ii) As Eq. Ijl6|) can be satisfied by /(r, p) — const, 
any solution of Eq. (|16|l can be shifted by a constant 
with no consequences for the physical quantities [at any 
measurement of the current one registers the difference 
between fluxes of incoming and outgoing particles and the 
isotropic part of /(r, <p) is cancelled out]. It follows from 
Eq. i|l(J|) that the distribution of particles impinging onto 
the diffusive region at the interface, Sn Vb ^ d (x = 0), is not 
affected by a global shift of the potential: Sn v ,,{x = 
0) = (dn° F /de)e [$(x = 0) - $(+oo)]. This is in lull ac- 
cord with the fact that one is free to shift the potential 
$ by a constant. 

In Sec. IIVI and [V] we choose $(+oo) = and corre- 
spondingly use f(x = 0, fb^d) = as the boundary con- 
dition for the function f(v,<p) at the interface between 
the diffusive and ballistic regions. 



IV. DIFFUSE EMISSION PROBLEM 



and the currentless "counterflow" g((,p), i.e., f(C,p) = 
/di-(C> p) + g(C, n )- At a distance about a few mean free 
paths I the counterflow g{C,p) approaches an isotropic 
distribution (generally, a nonzero constant). For cur- 
rentless counterflow particles injected into the diffusive 
region eventually return back to the interface. Hence the 
distribution function of particles emitted by the diffu- 
sive region, g{—p), should be determined completely by 
the distribution function of the injected particles, g{+p). 
Due to the linear nature of Eq. (|17|l . the g(±p) parts 
of the distribution function are related linearlj*ii by the 
angular Green's function S(p,p'): 



dn'S(fi,n')g(C,+u'). (25) 



Here d£l' stands for the angular integration in d dimen- 
sions, while the subscript + in the integral indicates 
that the integration is limited to the directions for which 
p! > 0. With the use of Eq. (|25|l one can reexpress the 
density p g (C) corresponding to <?(C, aO through the in- 
coming part of the distribution only: 



p q (() = / dnH(p)g(( 1+ p) 



where 



H{p) 



dn 



,S(fj,,fi') 



(26) 



(27) 



The Green's function S(p,p') satisfies a nonlinear in- 
tegral equation which can be derived as follows. Differ- 
entiate Eq. I|25|) with respect to £, express the derivatives 
dg/d( through the kinetic equation 1|17[) . and use Eq. 125|) 
to eliminate g(Ci~ u ) m favor of g^j+p). In result one 
gets 

J dQ> S(p,p>) + (28) 
= -^H(p) f dn'H{p')g{^+p'), 

id J+ 



In this section we present a solution of the classical 
problem of the diffuse emission (Milne problemii*i£) (see 
Fig.0J in the form convenient for the subsequent analysis 
of Sec. It is given for an arbitrary dimension d, but 
we are interested in the particular case of d = 2. 

Let us specify the notation for angle ip used through- 
out Sec. IIVI and M for each of the two regions, diffusive 
and ballistic, ip is chosen as an angle formed by the mo- 
mentum of an electron with a normal to the interface di- 
rected inwards to the corresponding region. Correspond- 
ingly, for p = cos ip we adopt the convention^ that for 
electrons propagating away from the interface p > (de- 
noted as +p in what follows), while p < (denoted — p) 
corresponds to electrons propagating towards the inter- 
face. 

The general solution of Eq. (|17|l can be written as a 
sum of the current carrying Drude flow [see Eq. (|18|l ] 




Diffusion 



Ballistics 



FIG. 4: Milne problem. The diffusive region is to the left of 
the interface, £ < 0. The ballistic region is to the right, £ > 0. 
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where Sd is a total solid angle in d-dimensions. As this 
relation holds for an arbitrary incident distribution the 
equation for the Green's function S(/J,, fJ) follows: 



1 

Sd~ 



H(jM)HQi'). 



(29) 



The counterflow does not carry current and therefore it 
satisfies 



dQ. \pg(+fx) - MS(-M)] = 0. 



(30) 



This implies that Green's function S(fi, fi!) should satisfy 
the condition 



(31) 



that can be verified with the use of Eq. l|29[l along with 
the relation J + dtt H(fi) = Sd- 

As it has been explained in Sec. II I II the geometry of 
the discussed setup is such that the distribution func- 
tion in the ballistic region does not contain a compo- 
nent describing particles propagating towards the in- 
terface, i.e., /(0, = 0. Therefore, in the solution 
/(C) A 4 ) = /Dr(C ; A*) + fl'Dr(C) A 4 ) one nas t° counterbalance 
the part / Dr(0, +fi) by a proper choice of <?Dr(0, = 
— /df(0, +/i). Then the emission outside the diffusive re- 
gion consists of the two contributions: /di-(0, — fi) from 
the Drude flow, and <7Dr(0, —fi) which originates from the 
compensating counterflow. As <?0 r (O, is known, the 
latter contribution can be found with the help of Eq. <I25|) . 
As a result, the emission in the Milne problem is given 
by the expression 



/ M (0,-/i) = J^-\M+- I dQ'fi'S^fx')}, 
Sd M J+ 



(32) 



which can be simplified using Eqs. (|29|l and (|31() together 
with the relation J + d£lfiH{fi) = Sd/Vd. Finally, the an- 
gular distribution of the particles emitted from the inter- 
face into the clean region is given by 



/m(0,-/x) = J^H{fx). 



(33) 



This result coincides with the one presented in Ref. [ill 
for the case d = 3. In Fig. [5] the function H(fj) and 
the angular dependence of the Drude flow distribution, 
~ cosy?, are plotted, both normalized to J ~ 1. As com- 
pared to the Drude flow, the diffuse emission distribution 
in the Milne problem has a qualitatively different behav- 
ior at large angles. Namely, the distribution flattens and 
a considerable part of the distribution is transferred to 
the large angles. 



f (£=0) 

0.4 
0.2 




7T/4 



ix/2 



FIG. 5: (Color online) The angular distribution, H(cosip), of 
the diffuse emission (red line) as compared to the Drude flow 
distribution (blue line) for J = 1. The function H(cosip) ^ 
at ip = 7r/2. 



V. DIFFUSE EMISSION IN THE PRESENCE 
OF SPIN-DEPENDENT SCATTERING AT THE 
N/SO INTERFACE. 

In this section we generalize the solution of the Milne 
problem to the case of spin-dependent reflection at the 
interface. We are mainly concerned with the influence of 
the strong reflection at the angles tangential to the in- 
terface on the distribution of the emitted electrons. We 
show that the effect of flattening of the angular distribu- 
tion is the robust property of the solution of the Milne 
problem which becomes even stronger in the presence of 
such a reflection. 

We now concentrate on the calculation of the diffuse 
emission in the presence of the spin-orbit interaction at 
the ballistic side of the junction. As compared to the 
consideration in the previous section the scattering at 
the interface now depends on the direction of spin of the 
incoming electron. On the other hand, electrons after 
a sequence of random scattering in the diffusive region 
return to the interface preserving their spin (see Fig. |SJ), 
as the impurity scattering of electrons in the diffusive 
region is assumed to be spin-independent. 




FIG. 6: (Color online) Milne problem in the presence of the 
spin-dependent scattering at the interface. The SO interac- 
tion in the diffusive region (£ < 0) is suppressed. 
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To describe such scattering, one should introduce a 
spin density matrix Q qq . In the diffusive region each 
component of the spin density matrix satisfies the kinetic 
equation l|17fl because in the normal region the Hamilto- 
nian is spin- independent. Consequently, one can apply 
the same angular Green's function S(n, //) as in Eq. (|25|l . 

We have to describe the counterflow at the interface in 
the presence of the reflection. For clarification purposes, 
let us imagine an auxiliary line, £ = —0, separating the 
diffusive region from the N-SO interface. (Actually, for a 
real device there can be a physical, not auxiliary, inter- 
val between the diffusion region and the beginning of the 
interface where a and B start to vary.) Similarly to the 
discussion of the Milne problem it is assumed that the 
electron distribution is given by a sum of the Drude flow 
and a counterflow. The incoming part of the counter- 
flow consists of the contributions: one, t>Dr(C = — 0, +/i), 
counterbalances the incoming part of the Drude flow, 
^Dr(C — ~0,+/i) = — f? Dr (C = — 0, while the other, 
v refl(( — — 0>+m)j i s determined by the flow incident 
onto the line ( = — 0. As it has been explained is Sec. IIIII 
and has been already used in Eq. we assume that 
there are no electrons incident onto the N-SO interface 
from the SO region. In the absence of electrons incident 
from the SO side of the interface the flow incident onto 
the line £ = — comes only as a result of the reflection 
at the N-SO interface [the latter circumstance explains 
our choice of a subscript for v re fi(( = — 0, +/i)]. There- 
fore the overall distribution of particles coming from the 
interface Q overall (( = -0, +fi) = v refl (C = -0, +fi). 

As a result of scattering inside the diffusive region 
fDr(C = ~~ 0> +/•*) transforms into a part of the outgoing 
counterflow t>Dr(C — ~ 0, — /i). Together with the Drude 
flow it yields the following contribution to the outgoing 
distribution g M (C = -0, -/i) = J{Vd/S d )H{p)a , where 
o o is the 2x2 unit matrix. The overall distribution of 
particles emitted from the line £ = —0 and incident onto 
the N-SO interface is given by 

BoveraiiiC = -0,-^) =J—H(n)a + - / dQ'S(fl,fl') 

A« J + 

xv refl (C = -0,+fi'). (34) 

The second term in this equation is generated by the in- 
cident part of the distribution, v r efl(C = ~0, +/i), which 
in turn is determined by the reflection of Q overa u (C — 
— 0, — A*) a t the N-SO interface. The relation between the 
incident and reflected parts of the distribution should be 
found from the solution of the scattering problem at the 
normal side of interface 

Vrefl(( = -0, +At) = K(p)Q overall {( = -0, -At)^(At), 

(35) 

where by 1Z we denote a 2 x 2 block of the scattering 
matrix corresponding to the reflection at the normal side 
of the interface. 

After substituting Eq. I|35|l in Eq. I|34|) one obtains a 



closed equation for Q overaH (n): 

1 If., 
Qoveraii(-V) =J-^-H(n)a + - / d&S(n,fi) 
V27T A* J+ 

xn(^)Q overaU {-^)n\^). (36) 

To analyze the effect of the reflection at the interface 
on the form of the diffuse emission, we study Eq. 13611 for 
the case of a smooth interface. In this limit the kinemat- 
ical aspect of the scattering at the interface is most pro- 
nounced and not masked by unnecessary complications. 
Namely, we assume that for electrons of the — chiral- 
ity the probability of the transmission is 1 at all angles, 
while for electrons of the + chirality the probability of 
the transmission is 1 for ip < ip c and otherwise; see 
Fig. |21 and the related discussion of the smooth interface. 

We will use a notation |x ± (±Ati*)) f° r spinors. The 
sign of the first argument indicates the sign of the pro- 
jection of the momentum of a scattering electron on the 
direction perpendicular to the interface, while i = ±1 
is the sign of the momentum component along the in- 
terface; the superscript ± denotes the chirality. Also in 
what follows the integration over tp will be performed as 

J dp = J2i=±i JLiX^/Vvl ~ M 2 )- I n this notation the 
reflection part of the scattering matrix for the smooth 
interface is given by 

K = 0, ip < <p c , 7 , 

= EJx + (+m,*)> (x + (-am)I, v>v c - {61) 

For each direction —At we introduce the orthogonal ba- 
sis n Q (— At, i) in such a way that n 1 is parallel to the in- 
cident momentum p; n 2 is perpendicular to the electron 
gas plane, n 2 || 1; and n 3 is directed along the vector 
of the polarization of the + chirality. For the analysis of 
Eq. (|36|) it will be convenient to use the four-component 
Bloch vector (s°,s) related to the matrix q as follows: 

s°(-fi,i) = -Tr<7°e(-AM)> 

s a (-fx,i) = ^Tv(n a a)g(-^i) 

q(-H,%) = S °a + ]T.s Q (-AM)(n a <T), (38) 

a 

where a is the unit matrix and <r are the Pauli matrices. 
In terms of the (s°,s) matrix Eq. I|36|l can be rewritten 
as a set of coupled equations 

1 1 pcosip" J / 

V2vr At y 1 Jo a/1 - At' 2 

4[s (-M',i) + s 3 (-M',i)] , (39) 

A* j Jo V 1 - M 

x(xW^I^-n"(-AM)|x+(+At',j)) 
x[ s °(V,j) + s 3 (-//,j)]. (40) 




FIG. 7: (Color online) Angles and directions of spinors in- 
volved in the derivation of Eq. I|41fl . For each angle if and if' 
there are two possible momentum directions corresponding to 
a different sign of the p z component. 

Remarkably, all components s a are determined by a sin- 
gle combination f£(-(Jt,j) = s°(-fj.,j) +s 3 (-fi,j). The 
combinations s°(—fi,j) ± s 3 (— describe the proba- 
bility for an incident electron to be in the ± chiral state. 
The representation of as s° ± s 3 is a consequence of 
the choice of the axis n 3 along the direction of the spin 
polarization of an electron in the state of the + chirality. 
The fact that the right-hand sides of Eqs. J55J| and l|4*U|) 
depend solely on f^(— //, j) has a simple physical rea- 
son: electrons of this chirality only are reflected at the 
interface. [In Eq. I|4(J|I the component s 1 is determined by 
A 4 ') i) without any feedback, and s 2 vanishes identi- 
cally in the chosen parametrization.] 

To proceed further we have to calculate 

(X + (+»',j)\k°- ■ n<*(-»,i)\ X + (+fi',j)) for both 
i, j = ±1. With the use of the relation that deter- 
mines the direction of the polarization of a spinor, 
<X + (+A*',j) klx + (+/Aj)) = n x +(+f*',j)> the discussed 
expression is reduced to ^ n x + (+n'-j) ' nQ ( — AM)- Then 
Eqs. Ij39|) and H40I ). can be rewritten in the form: 

1 1 fcosip* ii 

f+{-^ l = ±l)=J H^) + - / SM 

V27T /i Jo y/l - fi' 2 

xi{ [l-cos( p'±^)]/+(-//,+l) 
+ [l-cos(-^'±^)]/+(V,-l)}, 

(41) 

where angles y> and <p' are defined in Fig. [7\ For the 
difference 5f£(-(J,) = f^(-fi, -1) one gets 

a homogeneous equation 

^(-M) = - / . = , %,/)x (42) 

x \\- c °s(^' + y) + cos(^' - p)]Sf+ (-//). 

One can check that all eigenvalues of the right-hand side 
of this equation considered as a kernel of the transfor- 
mation are less than 1. Therefore we conclude that 
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tt/4 e c ji/2 



FIG. 8: (Color online) The displacement function /|o(#) of 
the " + " (shown in red) and "— " (shown in blue) polariza- 
tion components of the transmitted electrons for a = 0.1 as 
compared with the solution without reflection, 5 = (dashed 
line) . 

A*) = 0; i- e -; /n( — AM) is a symmetrical function 
with respect to i. Finally, this yields for Eq. 141fl 

V2tt A* Jo a/1 - m 2 

x (l-wO/^HA (43) 

Following the same route, it can be found from Eq. H4U|) 
the distribution of particles in the — chiral state: 

i i rcos<p c r / 

V^TT A* JO y/l — fl z 

x (1+MmO/n (-/"')■ (44) 

Equations (|43|l and (|44|l can be analyzed using the small- 
ness of cos(^ c , but they can be easily solved numerically. 
In what follows we will use the results of numerical anal- 
ysis for the distribution functions (— /z). 

Ultimately, we are interested in the distributions on the 
SO side of the interface. To achieve this goal, we have 
to connect the calculated distributions to the distri- 
butions on the SO side f§ . It follows straightforwardly 
from the Liouville's theorem that 

ifo^+Me) = /i(e,-AV), (45) 

where /z e and — //„ are two directions connected by the 
Snell's law, Eq. (g}. In the subsequent discussion the 
notations fg O (0) and f^(f) will be used instead of 

/sb( e > +A 1 e) an d /n ( e ' ~AV)' n °tice that we return to the 
definition of angles (f and 6* as it is given in Fig. ^ Fi- 
nally, in Fig. our result for /§o(0) is presented. 

We observe that the distribution of particles of the + 
chirality emitted from the diffusive region into the clean 
one is flat even in the presence of the spin-dependent 
reflection at the interface. Moreover, at large angles the 
function fs (9) is noticeably increased. 

As a result of the reflection at the interface there is 
a redistribution of the population of the particles of the 
+ and — chiralities. Consequently, the two chiral com- 
ponents of the current, j , change in the vicinity of the 
interface, and we are faced with the question of the con- 
servation of the total current j = j + + j~ . (We assume 
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for a moment the width of the diffusive stripe W to be 
unity, and until the end of this section will not distinguish 
between the density of current and the total current.) To 
check the current conservation we first calculated numer- 
ically the currents with the use of the distribution 
functions /^(e, — pb). We get that on the N side of the 
interface the current j + +j~ is equal to the total current 
j deep inside the diffusive region. Further on, it can be 
shown that the two parts of the physical current j do 
not change after passing the interface. To do this no- 
tice that in the SO region the chiral components of the 
physical current are 

3%o =Plo^(™s9f% (9)) (46) 

[for the definition of the physical currents see Eqs. i|19|) 
and H20fl ]. Having in mind Eq. (|45|l and the Snell's law 
in its differential form, pso cos9d9 — pn cos ipdip, the 
following relations for the ± chirality current components 
on the two sides of the interface can be obtained: 

J + = Pso ^2 [' % d0 cos 9f+ Q (9) = 
e 2 f p ° 

= Pn ~ n / dip cos ipf^(ip), (47) 
27r Jo 

and 

e 2 f eC 

3 = Pso 2^ J d9 cos 8 f SQ {9) = 

e 2 r /2 

= PNTT^T / dip COB (f fa (cp). (48) 
27r Jo 

For the + chirality component it has to be taken into 
consideration that the current of the particles impinging 
on the interface in the interval of angles ip c < ip < tt/2 is 
canceled out by the flow of the reflected particles. 

Thus we arrive at the conclusion that the two parts 
of the physical current do not change after passing 
the interface, i.e., = jg Q . Together with the fact of 
the current conservation in the N region this implies that 
the total current in the ballistic SO region is equal to the 
current deep inside the diffusive region. 

So far we limit ourselves to the diagonal elements of the 
density matrix. This was possible because in the N re- 
gion the density matrix is diagonal in momentum space, 
while the current operator is diagonal in spin space. The 
situation is different in the SO region where the current 
operator acquires the spin structure, see Sec. [nj In ad- 
dition, the density matrix becomes nondiagonal in mo- 
mentum space as an electron beam splits into two beams 
as a result of the refraction at the interface with the in- 
homogeneous spin-orbit interaction. We will not discuss 
this problem for the following reason. As it has been ex- 
plained in Sec. lIIII our analysis refers to the case when the 




spin filtered 
electrons 



FIG. 9: (Color online) A spin filter. Electrons emitted from 
the diffusive stripe at small angles to the interface are spin 
polarized and can be collected for subsequent spin manipula- 
tions. 

current is registered far away from the near-field zone of 
the orifice, i.e., at a distance L much exceeding the width 
of the diffusive stripe. At such distances the trajectories 
of electrons of different chiralities diverge from each other 
and therefore the nondiagonal in momentum space com- 
ponents of the density matrix become nonlocal in space. 
Since the current operator is local this nonlocality of the 
density matrix cannot show up in the current measure- 
ments. 



VI. SUMMARY 

We have analyzed the transport near the interface 
with the inhomogeneous spin-orbit interaction in terms of 
the spin density matrix. The present analysis has been 
performed under the assumption that the SO interac- 
tion in the diffusive region is absent (suppressed com- 
pletely). In fact, a much weaker condition is sufficient. 
It is enough that spin relaxation rate in the diffusive re- 
gion is controlled by the Dyakonov-Perel mechanism: 20 
1/t s ~ A| t < 1/r, where A 50 = 2a(x < 0)p F < 1/r 
is the spin splitting of the energy spectrum in the dif- 
fusive region, x < 0. Under this condition the spin re- 
laxation of the electrons during the propagation in the 
diffusive region after the reflection from the interface is 
negligible. (The other limit Ago ^ l/ r will be consid- 
ered elsewhere.) 

We have verified that the specific property of the so- 
lution of the Milne problem - the existence of the flat 
distribution with a large portion of electrons emitted at 
small angles to the interface - still holds in the presence of 
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the reflection at the interface. Moreover, at large angles 
the distribution function of electrons of the + chirality is 
noticeably increased, see Fig. |H1 As it has already been 
pointed out in Sec. [HI there exists an interval of the out- 
going angles, 9 C < 9 < ir/2, where only the + spin chiral 
component can penetrate. Together, these observations 
call upon to exploit the diffuse emission for the purposes 
of spintronics. 

In Fig. [5] a scheme of a device is presented which can 
operate as a spin filter for a current passing through the 
diffusive stripe confined between nonconducting areas A 
and B. Two additional barriers (or gates) are set at a 

distance L much larger than the width of the diffusive 

— ► < — 

stripe W. Within this geometry each collector, C and C, 
gets spin carriers emitted into the corresponding angular 
interval 59. When 59 < tt/2 — 9 C , particles of the + 
chirality only can get into the collectors. As a result, 
the currents inside each of the two collectors are spin- 
polarized, dominantly along the ± x directions. 

In the setup under discussion the orifice of the stripe 
acts as a source of a current with an angular distribution 
f% (9). The number of electrons of a certain chirality 
flowing in a direction 9 per unit time (i.e. the angular 
flux) is equal to 2^ (9)ded9, where the intensity 2^ (9) 
is related to f% (9) as 1% = j^yiPso f%o( 9 ) cos9W. 
The factor cos 9 W appears because we consider the flux 
of the particles outgoing from the orifice at an angle 9. 
The angular dependence of (9) is in full accord with 
the expressions for the density of the currents j^; see ex- 
pressions under the integrals in Eqs. (|47|l and (|48|l . As- 
suming the angular distribution of the emitted electrons 
to be practically flat, the spin-polarized current that can 



be collected by each of the collectors is ~ j 59 2 /8. As 
it has been explained in Sec. [H] the width of the angu- 
lar interval 59 = it/2 — 9 C can be as much as 2\fa. At 
this point the fraction of spin-polarized current reaches 
its optimal value ~ jot/2. All collected electrons have 
the same chirality that results in a very high level of spin 
polarization of the current in the collectors. A deviation 
from the perfect level of spin-polarization is only due to a 
small spread of the direction of motion of electrons within 
an angular interval 59. 

In Ref. a large spontaneous spin splitting has been 
detected in a gate controlled electron gas formed at a 
Ino.75Gao.25 AsZlno.75Alo.25As heterojunction. The re- 
ported splitting corresponds to a « 0.1. For such values 
of 5 one may expect a rather large angular interval 59 
that can be used for spin filtration; 59 w 36°. Under 
these conditions, the device proposed in Fig. [5] has the 
following specifications: a fraction up to 5% of the to- 
tal current is collected in C and is (almost) fully spin- 
polarized along the direction x, while the other fraction 
of 5% is collected in C and is spin-polarized along the 
direction — x. 

After filtration the spin-polarized current can be ma- 
nipulated similarly to the polarized light in optical de- 
vices. In particular, one can link the spin filter to the 
switch of the spin-polarized current discussed in Ref. Q]. 
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